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Abstract 

A material exhibiting a negative Poisson’s ratio is always one of the leading topics in materi¬ 
als science, which is due to the potential applications in those special areas such as defence and 
medicine. In this letter, we demonstrate a new material, few-layer orthorhombic arsenic, also pos¬ 
sesses the negative Poisson’s ratio. For monolayer arsenic, the negative Poisson’s ratio is predicted 
to be around -0.09, originated from the hinge-like structure within the single layer of arsenic. When 
the layer increases, the negative Poisson’s ratio becomes more negative and finally approaches the 
limit at four-layer, which is very close to the bulk’s value of -0.12. The underlying mechanism is 
proposed for this layer-dependent negative Poisson’s ratio, where the internal bond lengths as well 
as the normal Poisson’s ratio within layer play a key role. The study like ours sheds new light on 
the physics of negative Poisson’s ratio in those hinge-like nano-materials. 
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Virtually all known materials would undergo a lateral contraction when stretched lon¬ 
gitudinally and vice versa, which is protected by the conservation of volume under elastic 
loading. This gives a positive Poisson’s ratio (z/). In isotropic materials, the positive Pois- 


Q. 


son’s ratio is theoretically in the range from 0 (Cork) to 0.5 (Rubber) [If. If a lateral 
dimension expands during stretching or vice versa, the exhibited Poisson’s ratio is negative 
and the associated material is termed auxetic j^j]. Much intense interest in this counterin¬ 
tuitive feature stems from the pioneering discoveries that critical fluids, re-entrant polymer 
foams, colloidal crystals, laser-cooled crystals, unscreened metals, a-cristobalite, A7 struc¬ 
ture elements arsenic and bismuth, and laminates are found to be auxetic (sl-lt]] . Such an 
unusual mechanical property of u < 0 makes a material possible for applications in those 
special areas. For example, compression in one direction results in a shrink not an expansion 
in the transverse direction, demonstrating the mechanism of bulletproof vests in national 


defence [ 3 ]; The opposite situation—an expansion responds to a stretch—is the manifestation 
of artificial limbs in medicine js]. This sparks a surge in research activity in novel materials 
with negative Poisson’s ratios. 

In 2008, Hall and his colleagues found that negative Poisson’s ratios can appear in low¬ 
dimensional carbon sheets when multiwalled nanotubes are introduced 9]. This is a big step 
towards the design of nanostructured composites, artificial muscles, gaskets, and chemical 
and mechanical sensors. Naturally, it raises a question whether the negative Poisson’s ratio 
can emerge in other low-dimensional materials. In fact, it is difficult to measure this negative 
Poisson’s ratio in experiments as the observation of such a ratio is spurious jlo(]. But, one 
can always predict it theoretically [ll]. For example, monolayer of black phosphorus (BP) 
is reported to exhibit a negative Poisson’s ratio through using ab initio method 12]. This 
negative Poisson’s ratio is ascribed to the occurrence of hinge-like structure within a single 
layer of BP. Like BP, arsenic also possesses the hinge-like structures in the orthorhombic 
phase. More importantly, two groups independently demonstrate the thermal stability of 
its few-layer forms in this orthorhombic phase 14]. Their bandgaps are both layer- and 
strain-dependent and the gap values are around 1 eV. More surprising is that the carrier 
mobility can approach as high as several thousand cm 2 V -1 s -1 in these few-layer arsenic. All 
these make few-layer arsenic promising for future semiconducting applications. However, up 
to now, no study has been focused on the mechanical properties, in particular the negative 
Poisson’s ratio. Thus, at present, it is timely to check the possible existence of negative 
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Poisson’s ratios in few-layer orthorhombic arsenic. 

In this Letter, the negative Poisson’s ratio is reported to occur in few-layer arsenic through 
using first-principles calculations, which stems from the hinge-like structure within one sin¬ 
gle layer of arsenic. The negative Poisson’s ratio is about -0.09 at monolayer. When layer 
increases, the negative Poisson’s ratio becomes more negative and finally approaches the 
bulk’s value of -0.12 at four-layer. This layer-dependent negative Poisson’s ratio is demon¬ 
strated to closely follow the internal bond lengths, in particular the one perpendicular to the 
within layer direction, which indeed increase with layer. To better understand the underly¬ 
ing mechanism, a rigid mechanical model is proposed, shedding new light on the negative 
Poisson’s ratio in those hinge-like vdW nano-materials. 

First-principles calculations in this work are performed within the framework of density 
functional theory, as implemented in the SIESTA code 15]. We have used the generalized 


gradient approximation in the form of Perdew, Burke, and Ernzerhof functional 16]. The 
effect of van der Waals (vdW) interaction is taken into account by using the empirical cor¬ 


rection scheme proposed by Cooper [17|. Only the valence electrons are considered in the 
calculation, with the core being replaced by norm-conserving scalar relativistic pseudopo- 
18] factorized in the Kleinman-By lander form 19]. We have used a split-valence 


tentials 
doublc-C 
atoms 


}asis set including polarization orbitals with an energy shift of 100 meV for all 


20], The convergence is achieved when the difference of the total energies between 
two consecutive ionic steps is less than 10 ” 5 eV and the maximum force allowed on each 
atom is set to be 0.01 eV/A. 

We start our work from the orthorhombic bulk arsenic with space group Cmca, which 
is the same as BP. The conventional unit cell includes eight atoms, as shown in Fig. [Q]a). 
Each As atom within a single layer is covalently bonded with three As atoms, forming a 
puckered graphene-like hexagonal structure. This puckered structure is also called a hinge¬ 
like structure, which consists of two orthogonal hinges (atoms 456 and 612). This sets up 
the basis for those exotic properties in few-layer arsenic. The lattice constants are optimized 
to a = 4.70 A, b = 3.77 A, and c = 11.11 A, generating the internal parameters r ±2 = 2.58 A, 
r 34 = 2.56 A, 0i23 = 94.09°, and 6*234 = 99.05°. All these agree well with experimental and 
theoretical results jl3, 21]. A primitive unit cell and its Wigner-Seitz cell are illustrated in 
Fig. HKb). It is comprised of four inequivalent atoms, which is only one half of the bulk case. 
In principle, monolayer of arsenic (arsenene), as shown in Fig. [EI C )> can be obtained through 
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exfoliating its bulk counterpart. Once it is formed, some changes occur. The corresponding 
lattice constants a and b are changed to be 4.73 A and 3.71 A, respectively. In comparison 
with its bulk phase, a increases while b decreases. This gives a direct effect on the internal 
parameters: bond lengths r 12 and r ^4 are decreased to 2.53 A and 2.50 A, while the bond 
angle 0 12 3 is increased to 94.54° and leaves 6^34 = 100.36° nearly unchanged. As the thermal 
stability of a material is important for real device applications, here we confirm it from 
the cohesive energy. The calculated cohesive energy of orthorhombic bulk arsenic is 2.90 
eV/atom, which is smaller by 0.15 eV/atom than that of the A7 structure arsenic 22], This 
coincides with the fact that the A7 structure is the favored phase in experiment. However, 
when they are isolated into monolayers, they become equally stable as the cohesive energy 
difference is only ~0.01 eV/atom. This means monolayer of orthorhombic arsenic is possible 
in experiment. 

The calculated band structure of bulk phase is displayed in Fig. [1(d). It clearly shows 
that the valence band maximum (VBM) and the conduction band minimum (CBM) locate 
at the same crystal momentum Z point, demonstrating a direct bandgap semiconductor 
with nearly zero gap value. However, when it goes into monolayer, it behaves as an indirect 
bandgap semiconductor, as shown in Fig. [1(e). This is because that VBM appears along 
the T-X line and close to the X point, while CBM locates at the T point. The obtained 
bandgap is about 1 eV. The underlying mechanism of the bandgap transition from direct 
(bulk) to indirect (monolayer) is dominated by the mutual competition of the two interlayer 
bondings ri 2 and r 34 , see also the discussion in our previous work 13]. In the following, we 
pay attention on the mechanical properties of few-layer arsenic. 

When deformation is applied along the a; direction, the responding strain in the y direction 
is shown in Fig. [2(a). Note that the positive (or negative) e means a tensile (or compressive) 
strain. The obtained data (solid circles) behaves as a strongly nonlinear feature, which is 
well fitted by function of y = —v\x + z/ 2 x 2 + v^x 3 . The linear parameter iq is fitted to be 
0.35 and can be regarded as the linear Poisson’s ratio. Similarly, we can obtain the linear 
Poisson’s ratio in the z direction, as shown in Fig. [2](b) . The corresponding linear Poisson 
ratio is v = Vi = 0.13, which is nearly twice that in the y direction. This means that z 
direction is harder than y direction when they respond to the strain applied along the x 
direction. It should be noticed that in the hinge-like structure of arsenic (see Fig. Ha)), 
the deformation of z direction is dominated by the bond length r 34 and the bond angle #234 > 
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while that of y direction needs the changes of r\ 2 and # 123 . Perturbed by the same external 
field, such as the strain in our case, the bond angle is more easily affected than that of bond 
length. Thus, we can infer that the bond angle d 12 3 , in comparison with $ 234 , are largely 
changed when the deformation is applied along the x direction. This coincides with our 
first-principle calculations. 

If the deformation is applied along the y direction, the situation dramatically changes. 
In Fig. [2jc), e x linearly depends on e y . The fitted Poisson’s ratio is v = 1.07, which is 
much larger than the Poisson’s ratios when the deformation is applied along the x direction. 
This is the direct manifestation of anisotropic feature in mechanical properties, which is 
also the origin of anisotropic transport reported in our previous work 13]. Up to now, all 


the obtained Poisson’s ratios are normal and positive. However, when we check the strain 
along the z direction responding to the ^-direction loading, the law behaved as in Fig. [2]]d) 
is opposite. e 2 increases with e y , leading to an unusual Poisson’s ratio, namely v < 0. The 
fitted linear Poisson’s ratio is v = -0.093. Its magnitude is very huge if one remembers the 
negative Poisson’s ratio of around -0.027 in monolayer of BP 12]. This makes monolayer 
arsenic more suitable for applications in those special areas such as aerospace and defence 
where large negative Poisson’s ratios are requested. For monolayer arsenic, the internal bond 
lengths are larger than those of 2.42 A and 2.38 A in monolayer BP 12]. In principle, under 
the same loading, larger bonding would bear a stronger deformation. This explains why 
monolayer arsenic holds a larger (more) negative Poisson’s ratio, being consistent with that 


v mostly increases with atomic number Z 


3], 


It is well known that layer stacking is emerging as a new degree of freedom to tune the 
electronic properties of vdW hetero-materials [23]. However, its impact on the negative 
Poisson’s ratio is still missing. In the following, we will demonstrate how it affects the neg¬ 
ative Poisson’s ratios in few-layer arsenic. Once an additional layer is added, the negative 
Poisson’s ratio is changed accordingly as the effect of layer stacking is involved, as displayed 
in Fig. |3](a). In comparison with monolayer, the nonlinear feature of e z versus e y is signif¬ 
icantly enhanced for bilayer. This is represented by the larger fitted parameters of cubic 
terms, which are at least twice that of monolayer. In the vicinity of zero, L\ and L 2 are 
nearly identical. When the strain is beyond around ±4%, Li and L 2 largely separate from 
each other. This is because the two layers are no longer equivalent in bilayer due to the 
AB stacking. The separated feature in the stress-strain curves under larger strain reflects 
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the asymmetric structure of layers, confirming the validity of our simulations. The linear 
Poisson’s ratio is fitted to be about -0.14, which is increased by ~50% compared with that 
of monolayer. Based on the above discussion, the larger the bonding length, the larger the 
negative Poisson’s ratio. Compared to the bond length r 12 , the bond length r 34 is almost 
parallel to the z direction, which is directly related to the negative Poisson’s ratio. Due to 
the interlayer vdW interaction introduced by an additional layer, is increased to be 2.54 
A, resulting in a larger negative Poisson’s ratio in bilayer. 

In the case of trilayer, as shown in Fig. |3](b), the stress-strain curves L\ and L 3 are 
identical, leaving the curve L 2 deviated from L 3 and L 3 in the whole stress range. All these 
coincide with the symmetry of trilayer where layers L 3 and L 3 are symmetrically equivalent, 
but not for layer L 2 . This leads to two negative Poisson’s ratios. One is about -0.13 for the 
layers L\ and L 3 . The other one is about -0.098 for the layer L 2 , which is smaller than those 
of L\ and L 3 . This is because the length r 34 of layer L 2 is 2.51 A, which is smaller than 
the value of 2.55 A for layers L 3 and L 3 . When any more layer is added, the internal bond 
lengths in particular r 3 ^ will be not changed largely. For example, in the case of four-layer, 
the bond lengths r 3 4 are 2.55 A and 2.54 A for L\ (or L 4 ) and L 2 (or L 3 ), which are very 
close to the value of 2.55 A in the bulk arsenic. The obtained negative Poisson ratio of 
four-layer is about -0.121 (averaged value), as shown in Fig. |3](c). This value approaches 
the bulk’s value of -0.125, as shown in Fig. |3](d). This means that the quantum effect on 
the negative Poisson’s ratios is limited within four layers. 

To further understand the layer stacking on the negative Poisson’s ratios, we demonstrate 
it from the intrinsic puckered structure combined with the associated stress-strain curves 
e x versus e y . As pointed out above, the hinge-like structure is the origin of the negative 
Poisson’s ratios in few-layer arsenic. The reason is as follows. When few-layer arsenic is 
stretched in the y direction, the layer contracts in the x direction, protected by the normal 
Poisson’s ratios, as shown in Figs. |T](c) and IU(a)-|l](d). I 11 other words, atoms 3 and 4 as 
well as atoms 1 and 6 move inward along the x direction when stretched in the y direction 
(see Fig. Oja)). This directly causes the bond angles $234 and d 2 1 6 smaller compared to the 
initial values. Taking into account all the values of #234 and 0 2 1 6 (with/without strain) being 
larger than 90°, the layer thickness along the z direction (the projected distance of r 34 or 
r^) is increased, leading to a negative Poisson’s ratio. Therefore, the bond angles #234 and 
$216 being larger than 90° is another necessary condition for the negative Poisson’s ratios 
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in few-layer arsenic. This is also true for BP as the bond angle #234 is 97.64° in BP 12]. 
If a material possesses the bond angles 6*234 and 6*216 being smaller than 90°, the negative 
Poisson’s ratio will disappear. Future researches can test this prediction. 

Due to the layer stacking, the normal Poisson’s ratios are increased from 1.07 at monolayer 
(Fig. H](c)) to 1.19 at four-layer (Fig. 01(c)). This implies that the decrease of 6*234 or 6 * 2 i 6 
under the same stretching in the y direction is enhanced as layer increases. As a result, the 
negative Poisson’s ratio increases with layer. For four-layer arsenic, the normal Poisson’s 
ratio is very close to the bulk’s value (Figs. 0(c) and 01(d)). This is the reason why the 
negative Poisson’s ratio can approach the bulk’s limit at four-layer. This tells us that the 
effect of layer stacking goes into the negative Poisson’s ratio in the y direction which is 
indeed through the normal Poisson’s ratio in the x direction. 

In conclusion, the negative Poisson’s ratio is for the first time reported in few-layer 
arsenic through using first-principles calculations. The magnitude of negative Poisson’s 
ratio is about 0.1 at monolayer, which is about five times of magnitude larger than that of 
monolayer BP, suggesting extended applications in those special areas such as defence and 
medicine. When layer increases, the negative Poisson’s ratio become large (more negative). 
The limited value of around -0.12 is predicted at four-layer, which is very close to the bulk’s 
value. The underlying mechanism is demonstrated based on the layer-dependent internal 
bond length r 34 and the normal Poisson’s ratio of e x versus e y . The study like ours sheds new 
light on the layer-dependent negative Poisson’s ratio in those hinge-like vdW nano-materials, 
which will evolve into an active field. 
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FIG. 1: (color online), (a) Conventional unit cell of orthorhombic arsenic with lattice constants 
a , b, and c and internal atoms 1-6. (b) Primitive unit cell and its Wigner-Seitz cell (blue shaded 
configuration), (c) Top view of monolayer arsenic with rectangle showing the unit cell. Band 
structures of (d) bulk and (e) monolayer. 
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FIG. 2: (a) e y versus e x . The solid circles are simulated data and the line is fitted by function 
y = —u\x + U2X 2 + V3X 3 , with v\ = 0.35 as the linear Poisson ratio, = 1-00 and V3 = -1.00. (b) 
e z versus e x . The fitted linear Poisson ratio is obtained to v\ = 0.13, as well as = -0.48 and 1^3 
= 0.23. (c) e x versus £ y . Data are fitted to function y = — ux, with v = 1.07 as the linear Poisson 
ratio, (d) £ z versus £ y . Data are fitted to function y = — vix + V3X 3 , with the negative Poisson 
ratio v = v\ = -0.093. 
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FIG. 3: (color online), (a) e z versus e y for bilayer arsenic. Inset shows a four-layer arsenic with 
layers labeled as L\ — L 4 from bottom to top. The fitted functions for each layer are given as well. 
e z versus e y for (b) trilayer, (c) four-layer and (d) bulk arsenic. Other denotes are the same as in 
(a). 
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FIG. 4: e x versus e y for (a) bilayer, (b) trilayer, (c) four-layer and (d) bulk arsenic. The fitted 
functions are also given. 
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